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Many inflation models predict that primordial density perturbations have a nonzero three-point correlation 
function, or bispectrum in Fourier space. Of the several possibilities for this bispectrum, the most comm- 
mon is the local-model bispectrum, which can be described as a spatial modulation of the small-scale (large- 
wavenumber) power spectrum by long-wavelength density fluctuations. While the local model predicts this 
spatial modulation to be scale-independent, many variants have some scale-dependence. Here we note that this 
scale dependence can be probed with measurements of frequency-spectrum distortions in the cosmic microwave 
background (CMB), in particular highlighting Compton-y distortions. Dissipation of primordial perturbations 
with wavenumbers 50Mpc _1 < k < 10 4 Mpc _1 give rise to chemical-potential (p ) distortions, while those with 
wavenumbers 1 MpcT 1 < k < 50Mpc~* give rise to Compton-v distortions. With local-model non-Gaussianity, 
the distortions induced by this dissipation can be distinguished from those due to other sources via their cross¬ 
correlation with the CMB temperature T. We show that the relative strengths of the pT and yT correlations 
thus probe the scale-dependence of non-Gaussianity and estimate the magnitude of possible signals relative to 
sensitivities of future experiments. We discuss the complementarity of these measurements with other probes 
of squeezed-limit non-Gaussianity. 


The prevailing paradigm for the origin of the primordial 
density perturbations inferred from fluctuations in the cosmic 
microwave background (CMB) and from large-scale galaxy 
surveys is that they arose as quantum fluctuations during a 
quasi-de-Sitter phase (“inflation”) of expansion in the early 
Universe m. Inflation generally predicts a spectrum of pri¬ 
mordial perturbations that is nearly, but not precisely, scale- 
invariant, consistent with current measurements. It also gener- 
ically predicts that these perturbations should be very close to, 
but not precisely, Gaussian (i.e., have only small connected 
correlations beyond the two-point correlation function). The 
amplitude and detailed form of this non-Gaussianity varies 
considerably between models and is usually specified in terms 
of a three-point function, or “bispectrum” in Fourier space ED- 

The local model bispectrum, which peaks in the squeezed 
limit - £3 » k\ (where k\ < £2 < k^ are the three 
Fourier modes being correlated) appears in a number of sim¬ 
ple inflation models ||3] 0 and has thus become a standard 
workhorse. Local-model non-Gaussianity has been sought in 
the CMB 2H21, in large-scale structure (8], and through the 
scale-dependent biasing mm it induces in galaxy clustering 
at large scales. 

It has also recently been shown Ifl3l fl4) that local-model 
non-Gaussianity will induce spatial correlations between the 
CMB temperature fluctuation and chemical-potential (pi) dis¬ 
tortions as a function of position on the sky. These p distor¬ 
tions arise from dissipation of primordial perturbations with 
wavenumbers 50Mpc 1 < k < 10 4 Mpc at redshifts z — 
5 x 10 4 - 2 x 10 6 U3fl3. In the local model, the amplitude 
of these small-scale perturbations is modulated by the long- 
wavelength curvature perturbation that also gives rise to large- 
angle temperature fluctuations, thus inducing a pT correlation 
03131 . These pT correlations therefore probe local-model 


non-Gaussianity at wavenumbers 50 Mpc 1 < k < 10 4 Mpc~‘ 
far smaller than those accessible with CMB temperature fluc¬ 
tuations, galaxy surveys, or even future 21cm observations 

m. 

Still, the local-model bispectrum is just one of an infinitude 
of different types of bispectra. Large classes of bispectra have 
arisen from different ideas for inflation II19li22| , and no short¬ 
age of phenomenological parametrizations have been consid¬ 
ered. 

In this paper we point out that the correlation of CMB 
Compton-y distortions with the large-angle temperature fluc¬ 
tuation can be used, in tandem with the pT correlation, to 
probe different types of non-Gaussianity. Dissipation of pri¬ 
mordial perturbations with wavenumber^] 1 Mpc - 1 < k < 
50 Mpc 1 gives rise to Compton-y distortions to the CMB 
m. However, these spectral distortions have been largely 
disregarded because they are not easily distinguished from 
larger Compton-y signals from the intergalactic medium at 
low redshifts. Non-Gaussianity may, however, induce a yT 
correlation that could allow this dissipation-induced y distor¬ 
tion to be isolated from late-time effects (e.g., ||23l 1241 1 by 
means of its angular dependence. But even without being able 
to distinguish different contributions, one can still derive con¬ 
servative upper limits using yT correlations. Since y and p 
distortions probe different wavenumbers, the relative strength 
of the yT and pT distortions can be used to probe the scale 
dependence of primordial non-Gaussianity, extending previ¬ 
ous considerations of just the pT correlation (25). 


1 While modes with k < 1 Mpc~* also damp and create a y-distortion. after 

recombination the effective heating rate for these perturbations is much 

lower than at smaller scales mini. 
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Below, we first present a simple calculation, based upon 
the configuration-space description (rather than the Fourier- 
space description in prior work fi~3l 141) of local-model non- 
Gaussianity and the /jT correlation. This calculation illus¬ 
trates clearly that the /jT correlation arises from large-scale 
modulation of small-scale modes. The generalization to yT 
correlations is thus clear. We then parametrize the type of non- 
Gaussianity probed by the combination of \xT and yT correla¬ 
tions and estimate the sensitivity of future experiments. We 
close by discussing the complementarity of these measure¬ 
ments with other probes of squeezed-limit non-Gaussianity. 

To begin, the average jj and y distortions induced, with 
Gaussian initial conditions, by dissipation of primordial per¬ 
turbations are given by 


in terms of a Gaussian random variable r{x). We then write 
R(x) = < Rl(x) + R s (x), where Rl(x) is the part of R(x) that 
comes from long-wavelength Fourier modes and R s (x) that 
from short-wavelength Fourier modes, and similarly write 
r(x) = ri(x) + r s (x). By writing 

Rl+Rs = ?l + r s + ^/ni \ r L + 2r/,r s + r"j, (4) 

we infer that the small-scale curvature fluctuation in the pres¬ 
ence of some fixed long-wavelength curvature fluctuation 
Rl(x) is, to linear order in f n \. 


R s (x) = r s (x) 


1 + ^/ni^zX-T) 


(5) 


<F> ~ J' dlog kA R (k)W f j(k), 

(la) 

(y) ~ J' d log k A R (k)W y (k), 

(lb) 


where A} R (k) is the primordial curvature power spectrum, and 
Wfj(k) and W Y (k) define k-space window functions to account 
for the acoustic-damping and thermalization physics. Al¬ 
though more accurate expressions for the /.'-space window 
functions have been discussed Q7J |26l|27), here we will use 
the simple approximations, 

Wfj(k) * 2.3 [e -2 ** 7 *^ - e -2 * 2 ^], (2a) 

W y (k) * 0.4 [e -2 * 2 ^ - ^ 2k2/k ^y] , (2b) 

where k D (z ) - 4.1 x 10~ 6 (1 + z) 3/2 Mpc~' is the diffusion 
damping scale at redshift z. This quantity is evaluated at the 
initial and final redshifts at which /./ and y distortions are pro¬ 
duced, assuming that the transition between /r and y happens 
abruptly at z - 5 x 10 4 EH) [29]. Thus, k Djl = k n (2 x 10 6 ) =* 
1.1 x 10 4 Mpc -1 , k Djiy = k D (5 x 10 4 ) =* 46Mpc _1 , and 
&D,y = k n ( 1090) - 0.15 Mpc '. Evaluating the expressions 
in Eq. (JTJ) using A 2 (k) * 2.4 x 10~ 9 (k/ 0.002 Mpc' 1 )”'-', with 
n s = 0.96 (301, yields (/j.) * 1.9 x 10^ 8 and (y) * 4.2 x 10~ 9 , 
in good agreement with more detailed computations CD. 

Since the angle subtended by a causal region at the surface 
of last scatter is ~ 1°, we see in the CMB ~ 40,000 causally 
disconnected Universes. If primordial perturbations are Gaus¬ 
sian, then the amplitude A^(k sma n) of primordial perturbations 
on the small scales k sma n that induce spectral distortions will 
be the same everywhere. If, however, there is local-model 
non-Gaussianity, then the power spectrum A R (k s , x) for small- 
scale modes k s will differ from one causal patch centered at 
position x to another. The fluctuation will moreover be corre¬ 
lated with the long-wavelength curvature fluctuation R(x). 

This can be understood simply from the configuration- 
space description of local-model non-Gaussianity. The local- 
model curvature perturbation at position x is written, 

R(x) = r(x) + ^f n] r 2 (x), (3) 


Thus, the fractional change in small-scale power in a region 
of a fixed long-wavelength curvature fluctuation is 


d(R 2 ) 

~W) 


y/ni Rl&). 


( 6 ) 


We thus infer that, with local-model non-Gaussianity, the 
(fractional) chemical-potential fluctuation in a given region 
of the sky is given simply by the long-wavelength curvature 
perturbation in that region at the surface of last scatter. 

The same is true for the large-angle temperature 
fluctuation—it is determined primarily by the curvature fluc¬ 
tuation at the surface of last scatter and has magnitude 
AT/T ~ R/5. Therefore, for multipole moments £ < 100 that 
probe causally disconnected regions at the surface of last scat¬ 
ter, the cross-correlation between the (fractional) chemical- 
potential fluctuation A/r//r as diR 2 ^ / (R 2 ^ and the tempera¬ 
ture fluctuation AT/T has a power spectrum, 

C( T ~ 12/ nl Cj r . (7) 

This easily obtained result agrees with Ref. m, noting that 
their C*f is for the /./ fluctuation, rather than the fractional // 
fluctuation. The jj autocorrelation caused by non-Gaussianity 
is Cf 144 f 2 Cj T . Here we assumed that the trispectrum 
contributions are negligible. As our derivation clarifies, the 
f/T correlation arises from the squeezed limit of the bispec¬ 
trum, the part of the bispectrum that modulates small-scale 
power. For completeness, we include the large-angle temper¬ 
ature power spectrum, 


C 


TT 

t 


2n A r 
25 £({ + 1) 


6.0 x 10~ 10 
£(£ + 1) 


( 8 ) 


If the non-Gaussianity is scale-invariant, the yT and yy corre¬ 
lations are the same as the ///’ and //// correlations (all in terms 
of A y/y and A/r/yu) 

Cf * 12 UCf , Cf * 144 f 2 C[ r . (9) 


If, however, the non-Gaussianity is scale-dependent (see, 
e.g.. Refs. (3IH37)), then the value of f n \ that describes cor¬ 
relations between the long-wavelength (/q, < 0.01 Mpc 1 ) 
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modes responsible for large-angle CMB temperature fluc¬ 
tuations and short-wavelength modes (SOMpc 1 < k < 

10 4 Mpc ') responsible for the /j distortion, may differ from 
the value of f n \ that describes correlations between long- 
wavelength modes and 1 Mpc 1 < k < SOM pc 1 modes re¬ 
sponsible for the dissipation-induced y distortions. 

We therefore parametrize the scale-dependent non- 
Gaussianity in terms of a scale-dependent non-Gaussianity 
parameter f n \(k s ) defined by the squeezed-limit (k s = k 2 - 
k 3 » k\ = ki ) curvature bispectrum, 

B K (ki,k 2 ,k 3 ) ^ y/ni (k s )P K (k s )P K (k L ), (10) 

defined from 


(%,%%) = (2*) 3 <5 (3) (£i +k 2 + k 3 )B K (k, , k 2 , k 3 ) (11) 


for ki < 0.01 Mpc 1 . We then define two parameters, 
fl 1 = /m (k s - 7Mpc-‘) and f* = fjk s - 740MPW 1 ) 
to parametrize the non-Gaussianity on /j and y scales, re¬ 
spectively. Here, the f n \ parameter is evaluated roughly 
at the log-midpoints of the y- and /./-distortion k intervals, 
which defines the y- and /v-distortion pivot wavenumbers 
ky - 7 Mpc 1 and k fl - 740 Mpc ', respectively. Although 
the scale-dependence of f n \ is often modeled in the litera¬ 
ture as a power law in wavenumber k, this parametrization 
does not accommodate the possibility, which arises in curva- 
ton and multi-field models BSEBHSD, that / n i may change 
sign with scale. The above parametrization is therefore more 
general. We note that the k\ = k L <s k 2 — k 3 = k s 


bispectrum in Eq. (IOi is the squeezed limit of the com¬ 


mon bispectrum parametrization (e.g., ED), B<R(k\,k 2 ,k 3 ) - 
(6/5) \f n \(ki)P(k 2 )P(k 3 ) + perms], of scale-dependent non- 
Gaussianity. Our results are thus easily compared with many 
other results on scale-dependent non-Gaussianity. 

We estimate the detectability of the /jT signal following 
Ref. |fT3l assuming a noise power spectrum C^' n — w~ 
with W ; 1/2 « V47T (yt/min/ (//)) and ( max ^ 100. Note again that 
ours is a power spectrum for A/r/jU, as opposed to the power 
spectrum for /./ in Ref. oca, where jj mm is the smallest de¬ 
tectable jx monopole, which is estimated to be yi mm - 10 4 for 
PIXIE m and // min ^ 10~ 9 for PRISM |43]. We then obtain 
the smallest detectable (at ~ lcr) /jT correlation to arise for 


£“ 220 


( fl min \ 

HU 9 ] 


v> 


2x 10- 


( 12 ) 


Although the dissipation-induced average y distortion is ex¬ 
pected to be =: 5 times smaller than the fi distortion (see es¬ 
timate above), the experimental sensitivity to y is ^ 5-10 
times better l42l|44l . Assuming a PRISM detection limit of 
y m in — 2 x 10 l0 , we thus infer a comparable smallest non- 
Gaussianity parameter. 


£,“220 


/ finin \ 

1 (y) V 

\2 x 10 10 j 

\4x 10 9 / 


(13) 


that will engender a detectable yT correlation. These esti¬ 
mates are obtained assuming a roughly scale-invariant spec¬ 
trum of primordial perturbations. If for some reason the small- 
scale power-spectrum amplitude is increased in the y range 
(IMpW 1 < k < 50Mpc~') or /j range (50Mpc~' < k < 
10 4 Mpc~'), the smallest detectable and /^, respectively, 
will be decreased by a similar factor. In this case, the ho¬ 
mogeneous values of y and /j will also be increased Ifl7) . The 
detection limits obtained from the // and y autocorrelations are 
typically a few times weaker, since in contrast to the temper¬ 
ature cross-correlations noise dominates the ji distortion mea¬ 
surements. 

The /r and y spectral distortions in any given causal patch 
at the surface of last scatter come from the dissipation of 
small-scale modes within that patch. There is a finite num¬ 
ber N ~ (kA)' of such modes in this patch, where k is the 
largest relevant wavenumber GO 4 Mpc and 50Mpc 1 for ji 
and y distortions, respectively), and A ~ 100 Mpc the the size 
of the causal patch at the surface of last scatter. There will 
therefore be Poisson fluctuations of amplitude ~ /V 1/2 in the 
value of the // and y distortions in any such patch. The mean- 
square fractional /./ fluctuation in a patch of angular size 0 will 
thus be. 



m 

2n 


CAW A 2 ~ ec\ 


i~e-' 



where Wc is the window function for a circle on the sky of 
radius 6 , and the last (approximate) equality is the Poisson- 
fluctuation amplitude. We therefore infer fluctuations Cy ~ 
6 x 10 15 (^f/100) and Cf ~ 6 x 10~ 25 (^/100) well below 
those from measurement noise. 

Before the yT correlations considered here can be used to 
probe or constrain primordial non-Gaussianity, it will be nec¬ 
essary to consider the cross-correlation of the Compton-y dis¬ 
tortion from intergalactic gas in the late Universe with the 
contribution to the large-angle temperature fluctuation from 
the integrated Sachs-Wolfe effect. We anticipate that lensing 
reconstruction can be used to separate out this late-time yT 
correlation. We also anticipate that this contribution can be 
distinguished by a different ( dependence. 

An additional source of yT correlations could arise from 
the damping of primordial magnetic fields l48l . In spite of 
the large uncertainty in the amplitude of primordial magnetic 
fields, the scale-dependence is again generally expected to dif¬ 
fer from the one caused by non-Gaussianity. However, a more 
detailed study is required. 

The constraints to the squeezed-limit bispectrum from yiT 
and yT correlation considered here will be complemented on 
the longer- and shorter-wavelength ends by other measure¬ 
ments. Searches for non-Gaussianity in CMB fluctuations and 
in galaxy surveys probe wavenumbers k < 0.1 Mpc There 
are other searches for the scale-dependent bias that arises 
from non-Gaussianity SCUD- These probe squeezed-limit 
non-Gaussianity primarily on wavenumbers k < 1 Mpc -1 |[T2fl 
(those that are most important for determining the abundance 
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FIG. 1. The different k- space windows responsible for different ob¬ 
servables. The CMB anisotropies are visible at k < 0.1 Mpc -1 , where 
the smale-scale cut-off is introduced by Silk damping. Compton-y- 
distortions are formed by the damping of modes with 1 Mpc~' < 
k < 50Mpc~*, while yu-distortions are created by modes with 
50Mpc~* < k < 10 4 Mpc~'. Additional large-scale temperature 
power can be created by modes at 10 4 Mpc -1 £ k (45). Those with 
10 4 Mpc~* <, k <, 10 5 Mpc _1 can be constrained using the primor¬ 
dial helium abundance 14611471 . while modes with 10 5 Mpc -1 < k 
are erased before neutrino decoupling and thus would only lead ad¬ 
ditional large-scale temperature fluctuations. 


of the galaxies being correlated). Dissipation of acoustic 
modes with wavelengths 10 4 Mpc 1 < k < 10 s Mpc 1 pro¬ 
duce entropy in the primordial plasma after BBN [46, 471. 
The long-wavelength modulation induced by squeezed-limit 
non-Gaussianity on these scales will then give rise to a small 
isocurvature fluctuation correlated with the primordial adia¬ 
batic perturbation l46l . Finally, modes with 10 4 Mpc 1 < k 
can give rise to additional large-scale temperature fluctuations 
caused by non-Gaussianity (45]. Depending on the sign of the 
correlation, this could also produce a lack of power on large 
scales (which could also affect the values of cosmological pa¬ 
rameters inferred from these temperature fluctuations). 

We find it interesting that there are thus now prospects 
to probe the amplitude of squeezed-limit non-Gaussianity on 
a continuum of distance scales from the largest (~Gpc) ac¬ 
cessible to those on scales nearly eight orders of magnitude 
smaller, with CMB, large-scale structure, scale-dependent bi¬ 
asing, yT correlations, /jT correlations, and small-scale en¬ 
tropy production. This complement of measurements will 
thus allow the determination of the functional dependence of 
f n \(k s ), without necessarily assuming a specific parametriza- 
tion (e.g., power-law) for its scale-dependence. 
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